o 
o 

(N 



Dynamic bifurcations: 
hysteresis, scaling laws and feedback control 



OV N. Berglund 



School of Physics, Georgia Institute of Technology 
Atlanta GA 30332-0430, USA. 



O . 

<U , December 1, 1999 

Q 

Abstract 

^ ■ We review some properties of dynamical systems with slowly varying parameters, 

when a parameter is moved through a bifurcation point of the static system. Bifur- 
cations with a single zero eigenvalue may create hysteresis cycles, whose area scales 
in a nontrivial way with the adiabatic parameter. Hopf bifurcations lead to the de- 
layed appearance of oscillations. Feedback control theory motivates the study of a 
Q\ ' bifurcation with double zero eigenvalue, in which this delay is suppressed. 

On 

1 Introduction 

Many physical systems are described by ordinary differential equations (ODEs) of the form 

where x £ M n is a vector of dynamic variables and A £ l p a set of parameters. When 
modeling the system by the equation (|l]) we implicitly assume that the parameters are kept 
constant, while only the dynamic variables change in time. There are, however, situations 
in which the variables considered as parameters change slowly in time, for instance: 

• parameters which may be difficult to control, such as the temperature in a chemical 
reactor or the climate influencing an ecological system; 

• control parameters, such as the temperature difference in a convection experiment or 
the supply voltage of an electrical device, which are slowly varied in order to determine 
the bifurcation diagram experimentally; 

• slowly time-dependent forcings such as a periodic magnetic field acting on a magnet. 
It is thus important to understand the relation between solutions of the p-parameter family 
of autonomous ODEs (|l]), and the time-dependent system 



- = F(x,G(et)), (2) 

where the function G(t) : M. — * M p is given, e is a small parameter, and r = et is called 
slow time. It is convenient to rewrite (||) as the singularly perturbed system 

n t 

s—=F(x,G(r)) = f(x,T). (3) 
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There is a large literature on equations of this type, see for instance the textbooks jl], [|], 
the proceedings || and the paper Q for reviews of different approaches. 

We will focus here on relations between the solutions of the slowly time-dependent sys- 
tem (||) and the bifurcation diagram of the "frozen" system (|l]). It turns out that the dy- 
namics near nonbifurcating equilibrium branches is relatively simple and well understood. 
We summarize some of the relevant results in Section |2[ Bifurcations, however, lead to 
new interesting phenomena, some of which we will review in the next sections. 

• Bifurcations in which a single eigenvalue becomes equal to zero may lead to relaxation 
oscillations ||] , bifurcation delay and hysteresis || ^] . In Section [3|, we present some 
results on the e-dependence of the hysteresis area. 

• The Hopf bifurcation, in which a pair of complex conjugate eigenvalues crosses the 
imaginary axis, leads to the delayed appearance of large amplitude oscillations. In 
Section ^| we present some of Neishtadt's results ||, |9| on the determination of the 
bifurcation delay. 

Finally, in Section [5L we give a result on a bifurcation problem with a double zero 



eigenvalue, which arises in the context of feedback control pi, 11 1. 



2 Non-bifurcating equilibria 

If we naively replace e by in equation (||), we obtain the algebraic equation /(x,r) = 
which defines the equilibrium branches x = x*(r) of the frozen system. One can thus 
expect that for small positive e, some solutions of the slowly time-dependent system stay 
close to some equilibrium branch x*(t). This can be made precise in the following way. 
We consider the equation 

e^ = f(z,r) (4) 

under the following assumption. 

Assumption 1. There exist an interval Jgl, which need not be bounded, a neighbour- 
hood T> of the origin in M n , and positive constants ao, wo, d and M such that the following 
properties are satisfied uniformly in r G I. 

• The function f(x, r) : V x I -> R n is of class C k , k ^ 2. 

• There exists a function x*(t) : I — » T> such that /(a;*(r),r) = 0, with its derivative 
d T x*(r) bounded in norm by wo. 

• The matrix A(t) = d x f '{x* (r) , r) has eigenvalues aj(r), j = l,...,n, such that 
|Reaj(r)| ^ ao > for all j (hyperbolicity). 

• The function b(y, r) = f(x*(r) + y, r) — A{r)y is bounded by M||y|| 2 for ||y|| ^ d. 

Theorem 1. Under Assumption [1], there exist strictly positive constants Eq and c\, de- 
pending only on ao,u>o,d and M, such that if e sj eq: 

1. Equation (Q) admits a particular solution x(t) with 

\\x(t) -aj*(r)|| < cie Wei. (5) 

2. If k ^ 3, there exist functions Xj(r), j = 1, . . . , k — 2 and a constant Ck-i > such 
that 

k-2 

\\x(r) - [x*(t) +5> i z i (r)]|| < Cfc-xe*" 1 Vr G I. (6) 

3=1 
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3. If f is analytic in a complex neighbourhood of x*(t), there exist functions Xj(r), j 1, 
constants c,K>0 and an integer N(e) = 0(1/ e) such that 

N(e) 

\\x(t) - [x*( T ) + Y, ^Xj(t)}\\ < ce- 1/ft Vr G /. (7) 
i=l 

4- If f is periodic in r with period T, then x is also periodic with period T. 

5. If all eigenvalues of A(r) have a negative real part, there exist constants Mq,cq, k > 
such that any solution of (Q) with initial condition such that ||x(ro) — x*(ro)|| ^ cq at 
some tq £ I satisfies 

\\x(t) - x(t)\\ < M e- K ( r - r °)/ £ ||x(ro) - x(r )|| Vr G In [r ,oo). (8) 

6. If A has eigenvalues with both positive and negative real part, there exist local invariant 
manifolds on which the motion is either contracting or expanding. 

This result tells us that there exists indeed a particular solution tracking the equilib- 
rium x*(t) at a distance of order e. This particular solution is sometimes called adiabatic 
or slow solution. It admits asymptotic series in e which can be computed by substitution 
into (3). If the equilibrium is asymptotically stable, the adiabatic solution attracts nearby 
solutions exponentially fast. In this case, the relation 

lim x(t; e) = x*(t) Vr > r (9) 

implies that we may indeed take the formal limit e — > directly in (|j). 

Theorem [l| has a long history. In the asymptotically stable case, points 1. and 5. were 
originally proved in 13]. A different approach has been used in [14]. The exponential 



bounds in the analytic case are a result of an iterative scheme in H , an alternative approach 



can be found in [15]. The periodicity of solutions is a consequence of the implicit function 



theorem, and the computation of invariant manifolds in the hyperbolic case is explained 



in |0]. This result can be extended to periodic orbits |16 ]. 

Let us also point out that this result has an interesting consequence for the linear 
equation 

e^L = A(r)y, (10) 

which is similar to the Schrodinger equation as appearing in the adiabatic theorem in 
quantum mechanics 17, 18] . This equation also occurs if we linearize (||) around any par- 



ticular solution. The eigenvalues of the matrix A(t) G M nxn can be labeled by continuous 
functions ai(r), . . . , a n (r). 

Assumption 2. There exists a partition into two groups {a±, . . . , a p } and {a p +i, . . . , a n } 
such that the real gap 

7= inf |Re(aj(r) - aj(r)) I, (11) 



1 ^ i ^ P 
p+1 ^ j ^ n 



is strictly positive. 
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Theorem 2 (@). A ssume that A(t) is of class C 3 for r G I, and that Assumption [| 
holds. For sufficiently small e and t £ I, there exists an invertible matrix S(t,e) such 
that the change of variables y = S(t,e)z transforms (|To|) into 

ez = D{j,e)z, (12) 

where D(t, e) is block-diagonal, with one block of size p x p and eigenvalues aj (r) + 0{e) 
for j = 1, . . . ,p, and another block of size (n — p) x (n — p) and eigenvalues aj(r) + 0{e) 
for j = p + 1, . . . , n. The matrices S(t, e) and D(t, e) can be expanded into powers of e, 
up to exponential order in the analytic case. 

Since the connection between this result and Theorem |l| does not seem to be well 
known, we sketch the proof in Appendix 

Corollary 1. Assume that the eigenvalues of A(t,e) have uniformly disjoint real parts, 
that is, 

inf |Re(ai(r) - a^r)) I > 0. (13) 

1 ^ i<j ^ n 

Then equation (|10|) can be diagonalized by a change of variables y = S(t, e)z, and thus its 
principal solution can be written in the form 



U(t,t ) = S(t,s) 



/ e ai(r,T )/£ 

V o 




S(T ,e)-\ (14) 



where 



<Xj{T,T ) = I aj ( S )ds + 0(£). (15) 

J TO 

3 Single zero eigenvalue: hysteresis and scaling laws 

Theorem |] breaks down if some eigenvalues of the linearization A(t) = d x f(x*(T),T) 
cross the imaginary axis, that is, in the case of a bifurcation. In this case, new phenomena 
may occur. As noted in the center manifold theorem can be applied to the equa- 
tion ^-(x,t, e) = (f(x,r),£, 0) to reduce the dimension of the system to the number of 
bifurcating eigenvalues. 

The simplest case occurs when a single eigenvalue vanishes at the bifurcation point. 
We write the reduced equation in the form 

dx 

e ^ = /( 2; ' r )' ^K. (16) 

We shall assume that a stable equilibrium branch x*(t), existing for negative r, bifurcates 
at the origin. This requires that /(0, 0) = d x f(0, 0) = 0. 
Let us consider a few examples: 

1. In the saddle-node bifurcation, the stable branch meets an unstable branch, and no 
equilibrium exists for positive r. This results in the adiabatic solution jumping to some 
other attr actor. Jump phenomena were first studied in |]|, and later in (2^, § 0. 
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Figure 1. Comparison between solutions of equations with slowly moving parameters 
(light curves) and the bifurcation diagram of the associated frozen system (heavy curves, 
where full curves represent stable equilibria, and broken curves unstable ones), (a) The 
bifurcation diagram of equation (|l7]) contains two saddle-node bifurcations. Solutions 
converge to a hysteresis cycle with area satisfying A(e) — -4(0) ~ e 2 / 3 . (b) The bifurcation 
diagram of equation ( [L8| ) contains two pitchfork bifurcations. Solutions track the upper 
branch for increasing A and the lower one for decreasing A, creating a hysteresis cycle with 
area satisfying -4(e) — -4(0) ~ e 3 / 4 . 



2. In the transcritical bifurcation, a stable and an unstable branch exchange stability. It 
has been shown in plj that generically, solutions will track the outgoing stable branch, 
although it occasionally happens that the unstable branch is tracked for some time 



m 



3. In the pitchfork bifurcation, the stable branch becomes unstable, and two new stable 



branches are created. It has been shown in [23] that generically, solutions will track 
one of the two outgoing stable branches, where the choice of the branch depends on the 
geometry. If the original equilibrium does not depend on r, however, the bifurcation 
may be delayed || 24]. 

These local properties of bifurcations have interesting global consequences if the pa- 
rameter is varied periodically. The bifurcation points may be arranged in such a way that 
different branches are tracked for increasing and decreasing parameter, and thus there is 
hysteresis (and relation @ no longer holds). A classic example is the equation 



dx 
"dr" 



ar + A, 



A 



SUIT, 



(17) 



which leads to a familiar hysteresis curve (Fig. 0a). The curves x(t) display so-called 
relaxation oscillations in which slow and fast motions alternate. Another example is 



dx 
: dr" 



(x + A)(l - A 2 - (x + A) 2 ), A = 2sinr, 



(18) 



see Fig. ||b. In combination with other global mechanisms, such as inertia in the equation 

d 2 x 



dr 2 



-7£— + \(t)x - x , 
dr 



A 



SUIT, 



(19) 



hysteresis may also lead to a chaotic motion || 25 1. 

Another interesting effect of bifurcations is that equilibrium branches are no longer 
tracked at a distance of order e, but at a distance scaling with some fractional power of 



s. This results in hysteresis areas scaling in nontrivial ways with e [26, 27, 28 1. Similar 
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scaling laws have been observed in ferromagnets [29, 3C, and references therein], and they 
continue to be a domain of active research. 

We present here a method to determine scaling laws in an easy, geometrical way || f|j . 
The key idea is to understand the behaviour of the variable y = x — x*(t), which satisfies 
an equation of the form 

ep- = a(T)y + b(y,T) + £ W (T), (20) 
dr 

where a(r) = d x f(x*(r),T), w(t) = — d T x*(r) and b(y,r) = 0(y 2 ). 

Notation 1. We write x(t, e) ~ y(T,e) if there exist constants c± > 0, uniform in r and 
e, such that c_y(T, e) ^ x(t, e) ^ c + y(r, e) for small r and e. 

Assumption 3. 

• The function /(x,r) is of class C k , k ^ 3, for ^ x ^ d and tq ^ t < 0, and admits 
the Taylor expansion 

f(x,r)= Yl CnmX n T m + ^ R nm (x, r)x n T m , Coo = Clo = 0. (21) 

n+m<k n+m=k 

• For tq ^ r < 0, there exists a function x*(r) ~ |t| 9 (q > 0) such that /(x*(r),r) = 
and f(x, t) < for x*(r) < a; < d. 

• The function a(r) = d x f(x*(r),T) satisfies ci(t) ~ — |r| p , where 

p = min \q(n — 1) + m \ (n + m < k and c nm ^ 0) or n + m= fc}. (22) 
n ^ l,m ^ 

The numbers q and p can be determined geometrically by Newton's polygon, which is 
defined as the convex envelope of all points in (n, m)-plane such that either n + m < k 
and c nm ^ 0, or n + m = k. Then — q is the slope of a segment of the polygon, and p is 
the ordinate at 1 of this segment. The condition (]22|) is generically satisfied, as can be 



seen by inserting x*{r) ~ \t\ q into the Taylor expansion (21). 

Theorem 3 (@). Under Assumption^, any solution starting sufficiently close to x*(tq) 
reaches the O (e) -neighbourhood of x*{t) at a time t\ = tq + C(e|lne|) and satisfies 



x(t) - x*(t) 



, , „ forn^T^ - E V(P+1) 

| r |P+i-5 (23) 

£ q/(p+l) f Qr _ £ 1/(P+1) <; T <; 0. 



The area between x(r) and x*(t) scales as 

• e fo+l)/(p+i) if q< p 7 



e|lne| ifq = P, (24) 

e ifq>P- 



We sketch the proof in Appendix ^. It is more difficult to study the behaviour of 
solutions after the bifurcation, where there is no such general result. One can, however, 
obtain much information by studying the dynamics of the distance between x and the 
various equilibrium branches emerging from the bifurcation point, using similar methods. 

For the branch of the pitchfork bifurcation with vertical tangent, we have q = | and 
p = 1, which implies that this branch is tracked at a distance of order e 1//4 , and the area 
scales as e 3 / 4 . For the saddle- node bifurcation, we have q = p = \. The branch is tracked 
at a distance of order e 1 / 3 , and the jump turns out to be delayed by a time scaling as 
£ l/(p+l) = £ 2 / 3 5 which accounts for the scaling law in Fig. |l|a. 
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4 Hopf bifurcation and bifurcation delay 



We consider now the equation (|j) in the case where a pair of complex conjugate eigenvalues 
crosses the imaginary axis from left to right. The frozen system generically displays one 
of the following two behaviours: 

• in the supercritical Hopf bifurcation, a stable periodic orbit is created, 

• in the subcritical Hopf bifurcation, an unstable periodic orbit is destroyed. 

When the parameter is slowly varied through the bifurcation point, one would expect 
the appearance of oscillations with a slowly increasing amplitude in the first case, and a 
jump transition in the second case. This, however, does not occur if the system is smooth 
enough: the moment at which the trajectory departs from the unstable equilibrium turns 
out to be delayed with respect to the moment of the bifurcation (Fig. la). This bifurcation 



delay was first described in [31], analysed rigorously by Neishtadt |, ||, and observed in 
various physical systems [32, |33|]. We now state a version of Neishtadt's result on the 
computation of the delay. 

Assumption 4. There exist an open interval I £ K containing and a neighbourhood 
T> of the origin in R n such that 

• The function f{x, r) is analytic in a complex neighbourhood of / x T>. 

• There is a curve x*(t) : I — > M n with /(x*(t),t) = 0. 

• The matrix A(t) = d x f(x*(T),T) has two eigenvalues a(r) ± iu;(r), where a(0) = 0, 
a'(0) > 0, w(0) / 0, and a(r) has the same sign as r; all other eigenvalues of A(t) 
have a strictly negative real part. 

For r in a complex neighbourhood of /, we can define the function 

$( r ) = f [a(s)+iu(s)]ds. (25) 
J o 

For negative to G /, we define the function 

n(r ) = sup{r | Re*(s) < Reip(T ),T < s < r}. (26) 

T>TQ 

We have II(to) > for tq < because a(s) is negative for negative r, which implies 
that Re x I / (r) is decreasing for negative r. Moreover, we have lim T0 ^o- n( r o) = and 
lim To ^o- n'(ro) = —1. For instance, if a(r) = r, then II(to) = —tq. 



Theorem 4 (|8[). For any strictly negative tq G / there exist constants £o> c O) c i > and 
a continuous function r(e) with lim e ^o r ( e ) = such that for e < eq, any solution with 
initial condition such that ||x(to) — x*(to)|| ^ cq satisfies 

\\x(t) - x*(t)\\ ^ ci£ for t + r(e) ^ r ^ min{r + , n(r )} - r(e). (27) 

The so-called buffer time t + > depends only on the behaviour of the eigenvalue a{r) + 
iw(r) in the complex plane. 

The computation of the buffer time r + is explained in ||, see also ||. Roughly speak- 
ing, it is the largest real time which can be connected to the negative real axis by a path 
with constant Re\E'. It is usually of the order of lo(0). The remarkable fact is that unlike 
for the pitchfork bifurcation, the bifurcation delay exists also when the equilibrium x*(t) 
depends on r. 
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FIGURE 2. (a) The slow passage through a Hopf bifurcation (occurring at the origin) leads 
to the delayed appearance of oscillations. The delay depends on the initial condition, but 
admits a saturation value, called buffer time, (b) The delay is suppressed by an appropriate 
feedback control, which creates a bifurcation with double zero eigenvalue. Solutions track 
a stable equilibrium branch emerging from the bifurcation point. 



5 Double zero eigenvalue and feedback control 

Jumps and large amplitude oscillations may have catastrophic consequences for a device 
undergoing bifurcations. This has lead some researchers to design feedback controlled 
systems. A simple affine control system takes the form 

dx 

— = f(x,X)+bu(x,\), (28) 

where f(x,X) is the uncontrolled vector field satisfying Assumption |, 6 G M" is a given 
vector, which defines the direction in which the system can be steered, and u(x, A) £ M 



is the scalar feedback control. In previous works [34, |36|, feedback controls have been 
designed in order to make the bifurcation supercritical, so that orbits departing from the 
unstable equilibrium have a bounded amplitude. 

The results of the previous section show however that in the case of a dynamic Hopf 
bifurcation, a jump transition occurs even in the supercritical case. To avoid such a 
jump, one may try to suppress the bifurcation delay by moving the imaginary parts of the 
bifurcating eigenvalues to zero, in order to push the buffer time to zero. This will produce 
a bifurcation with double zero eigenvalue. Moreover, the control should ensure that the 
system is nonlinear ly stable at the bifurcation point, in such a way that stable equilibria 
are created which attract the solutions. 



In 1 1C ] we construct a feedback control which satisfies these requirements. We then 



show that the dynamics can be reduced to the two-dimensional equation 

£ t = V (29) 
= h(t)x + 2a{r)y + ^{t)x 2 + 5{r)xy - x 2 y - x 3 + 0({x 2 + y 2 ) 2 ) + eR(x, y, r, e), 

dr 

where a(r) is the same as in Assumption the functions ^(r),7(r) and 5(t) vanish 
at r = and i?(0,0, r, 0) is directly related to the drift of x*(t). The linearization of 
(]29| ) at the origin has eigenvalues a(r) ± /i(r) 1//2 . The frozen version of this system is 
a codimension-four unfolding of the singular vector field (y, —x 2 y — x 3 ) which has been 



studied in detail 37] 
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Theorem 5 ( fllOfl ). Assume that the right- hand- side of (|29| ) is analytic, a'(O) > and 
//(O) > 0. There exist positive constants d, T, M, k and a neighbourhood J\f of the origin 
in 1R 2 with the following property. For every To £ [— T, 0), there is a constant c\ > such 
that for sufficiently small e, any solution of (29) with initial condition (x,y)(ro) G N 
satisfies 

\x{t)\ Me\r\-\ \y(r)\ ^ Me\r\~ l l 2 , n(e) ^ r ^ - (e/d) 2 ' 3 , (30) 
\x{t)\ ^ Me 1 / 3 , \y(r)\ ^ Me 2/3 , ~{e/d) 2/3 {e/d) 2 / 3 , (31) 

where t\(e) = to + cie|lne|. //, moreover, the relations /x'(0) > 2a'(0) and R(0, 0, 0, 0) 7^ 
/io/d, i/ien for (e/d) 2 / 3 ^ r ^ T we have 



\x(t) - x+(t)\ ^ M^t- 1 + e 1 / 2 T~ 1 / 4 e~ KT / £ ], 
\y(r)\ ^Mfer-^+^/V^e-^], 

where 



(32) 



x (r) = J^ + °( r )' ##(0,0,0,0) >0, 
\-VA* + 0(r), i/ #(0,0, 0,0) <0 

are equilibria of (|29|), i.e., £/ie right-hand side of (^) vanishes when x = x+,y = and 
e = 0. 

This behaviour is illustrated in Fig. |2|b. The important fact is that instead of tracking 
the unstable equilibrium (represented by the origin) for some time, the trajectory smoothly 
tracks the equilibrium x + (r), a feature that can be used to detect the fact that a bifurcation 
has occurred. The behaviour is more complicated if /u'(0) < 2a'(0) Jl0|, O], 
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A Proof of Theorem |2] 

If A(t) is of class C k and satisfies Assumption 0, there exists a matrix Sq(t) of class C k 
such that Sq 1 ASq = Dq is block-diagonal p9f] . The change of variables y = So(T)yo 
yields the equation 

e^ = M r)y , M r-,,) = - sS^ = . (34) 



This equation can be transformed into the block-diagonal form ( |12[ ) by the change of 
variables yo = S\(t)z if we impose that S\ satisfy the differential equation 

e^=A (r)S 1 -S 1 D(r). (35) 
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We will look for a solution with matrices of the form 



l p eS 12 \ fDx 



CJ- D(T) = W D2 J- < 36 > 

where 5 i2 G R px(n - !,) and S al G r("-p) x p. Substitution into © leads to 
diS' 

e — : — = A\2 + ^4.iiiSi2 — 5"i2^422 — £ 2 S12A21S12, D2 = j4 2 2 + £ 2j 42i»S'i2 (37) 
dr 

and similar relations for S%i,Di. Let L be the linear operator L : X ^ A\\X — XA22- 
For e = 0, the eigenvalues of L are exactly given by Oj — dj, 1 ^ i p, p + 1 ^ j ^ n 
[38, |39[ |. Assumption [2] implies that L is invertible (its inverse can be represented as a 



complex integral |4(|]), and by the implicit function theorem, the right-hand-side of 
vanishes for S12 = S*(t) = —L~ 1 A\2 + 0(e). The Frechet derivative of ( |37|) around S* is 
equal to L + C(e 2 ), and thus, by Assumption ||, S* is hyperbolic for small e. Theorem [j] 
can then be applied to show the existence of a solution S\2{r) = S*(t) + 0{e). □ 



B Proof of Theorem |3| 

The assertion on t\ is a consequence of Theorem |]. Let \i = ^-j- and v = ^j-j-. We first 
consider equation ( |20[) for n ^ r ^ — e v . By assumption, we have a(r) ^ — a+|r| p for 
some a + > 0. Let bo > 0. Since y(ri) = 0(s), there exists by continuity a time r* G (ti, 0] 
such that y(r) ^ a+|T| 9 /26o for ri ^ r ^ r*. Moreover, we may assume that either 
r* = —e u or y(r*) = a+\r\ q /2bQ. Using the Tayor expansion of /, one can show that for 
r G [ri,r*], b(y,r) ^ &oy 2 | r l p ^ <? ] provided 60 is large enough. Then for r% ^ r ^ r* 

-^M^ + ^M 9-1 - (38) 

dT 26o 

This linear equation can be solved explicitly, to prove that y(r) ^ c+e|r| 9_p_1 on the 
time interval under consideration, for some c + > 0. Let us write r* = — e u ~ s ,5 ^ 0. 
Then y(r*) < c+e 5 ^ 1 ) It*) 9 . By definition of r*, we have either r* = — e", or y(j*) = 
a+\T*\ q /2b , and then e 5 ^ 1 ) ^ a+/26 c+, which shows that r* ^ - (26 c+/a+) I/ ^ « -e v . 
One obtains a lower bound in the same way, which proves that y{r) ~ e|r| 9 ~ p ~ 1 for 

For t* ^ r 0, the rescaling y = e^z, r = e^cr gives 

dz 

— = a(a)z + b(z,a)+w{a), (39) 
der 

where a(a) = e v ~ x a{e u a) » — |<r| p , u)(cr) = e u ~ li w{e v a) M 9-1 , and one can show with 
the Taylor expansion that b(a,z) = e v ~ ii b{e il z^£ v o) is smaller than a constant which 
does not depend on e. Using the fact that d a z ^ —a_\a\ p z — b z 2 , we obtain by solving this 
Bernoulli equation that z is bounded below by a positive constant uniform in e. Moreover, 
the hypothesis f(x*(r) + y,r) < for y > shows that 5(<x)£ + b(z,a) < for z > 0, 
which yields the upper bound d CT z ^ u;q|<t| 9 ~ . This shows that z(a) is also bounded from 
above by a constant independent of e. This proves that y(r) as e M for r* ^ r ^ 0. □ 
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